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Quantum properties of the probes used to estimate a classical parameter can be used to 
attain accuracies that beat the standard quantum limit. When qubits are used to construct 
\ a quantum probe, it is known that initializing n qubits in an entangled "cat state," rather 

than in a separable state, can improve the measurement uncertainty by a factor of 1/ \fn. 
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^vq ' We investigate how the measurement uncertainty is affected when the individual qubits in 

a probe are subjected to decoherence. In the face of such decoherence, we regard the rate 
] ^ ■ R at which qubits can be generated and the total duration t of a measurement as fixed 

resources, and we determine the optimal use of entanglement among the qubits and the 
. resulting optimal measurement uncertainty as functions of R and r. 
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2 ' I. INTRODUCTION 

cr 

This paper considers the question of quantum limits on estimating the value of a parameter 
that influences the state of a physical system . We call this system, which is 

intrinsically quantum mechanical, a probe, because it is used to probe the value of the parameter. 



The accuracy with which the parameter can be estimated is determined by the initial quantum 
state of the probe, the type of interaction by which the parameter influences the state of the probe, 
and the readout measurement that is used to extract information from the probe. 

In this paper the probe is always a collection of n qubits. We assume that there is sufficient 
control over the probe qubits to initialize the probe in any separable or entangled pure state. We 
^ ■ let g denote the parameter we are estimating. The effect of g on the jth probe qubit is described 
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by the Hamiltonian 



Hj(g) = \go- z -,j ■ (l.l) 



Here a z -j denotes the Pauli z operator for the jth qubit; similarly, a x -j and a y j denote the other two 
Pauli operators. The parameter g, which has units of frequency (we choose fi = 1), is a coupling 
strength. The Hamiltonian (jl.ip generates a rotation about the z axis of the qubit's Bloch sphere. 
The overall influence of the parameter on the probe is given by the Hamiltonian 

n 1 n 

H = y £H j (g) = ~g y E* t . J . (1.2) 

3=1 3=1 

The value of g is to be deduced from the change in the state of the probe. For simplicity, we assume 
that the probe qubits do not have any free Hamiltonian evolution. 
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Giovannetti, Lloyd, and Maccone jij] have analyzed a general scheme of this type. Their theoret- 
ical framework involves estimating a dimensionless parameter <p, introduced on the probe through 
a unitary transformation U = e~ lhip , with h = Y^]=\ hj- They do not identify h as a Hamiltonian; 
rather it is treated as an arbitrary operator that is the generator of translations in the variable cp. 
The connection between our scheme and theirs is established by identifying h = H/g (hj = Hj/g) 
and <p = gT, where T is the time for which the Hamiltonian (|1.2[) acts on the probe. 

The chief objective of our analysis, following [hJ, is to expand the discussion in [9] by in- 
vestigating how decoherence impacts the accuracy with which the parameter can be determined. 
Thus we assume that in addition to the Hamiltonians Hj(g), the probe qubits are subject to other 
influences that can lead to decoherence. For the decoherence models we consider, the effects of 
decoherence manifest themselves at a readily identifiable rate, which we denote as 7. To make the 
analysis meaningful, we must impose additional constraints on the probes, since we can always 
make decoherence irrelevant by estimating g using a procedure that is completed in a time much 
shorter than the time 7 -1 over which decoherence has a significant effect. Thus we assume that 
qubits are made available and initialized into probes at a rate R. What we have in mind is that 
each probe is assembled in a time n/R and is then sent immediately through a quantum channel, 
where it is subjected to the Hamiltonian fll.2|) for a time T. If we use v probes to estimate g, so 
that the total number of qubits is 

N = vn , (1.3) 

the total time required is 

r = vn/R + T , (1.4) 

provided that the quantum channel can accommodate more than one probe at a time. In our 
analysis, we assume that the parameter must be determined in the fixed time r (r _1 can be 
thought of roughly as the bandwidth over which a time- varying g is estimated), that the qubit 
supply rate R is a fixed resource, and that the decoherence rate 7 is a constant. We vary the 
interaction time T and the number of qubits in each probe, n = R(t — T)/v, to achieve the best 
accuracy in determining g. 

A measurement scheme of the sort discussed here appears in slightly modified forms in several 
problem of practicaf importance, such as Cock ^chroni^ion Q Q Q Q, reference- 



frame alignment [ttI. [l8 1 , p hase estimation 0, 0, 21], frequency measurements 10, 22, 23], and 
position measurements 



The accuracy with which g can be estimated is closely connected to the distinguishability 
of neighboring states of the quantum probe. This connection is quantified by the generalized 
uncertainty relations formulated by Braunstein, Caves, and Milburn As in 0], we use these 
generalized uncertainty relations to describe the optimal accuracy of parameter estimation. 

In Sec. HH we review the formalism of generalized uncertainty relations in the forms suitable 
for our analysis and discuss briefly aspects of the assumptions we make about resources and time 
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scales in our measurement protocol. Section IIIII reviews the accuracy that can be achieved in the 
absence of decoherence, and Sec. HV] investigates how the achievable accuracies are affected by a 
general qubit decoherence process. The final section provides a short discussion of our results. 



II. GENERALIZED UNCERTAINTY RELATIONS 



As a consequence of spending a time T in the quantum channel, the quantum state of the probe 
changes from an initial state po to a final state 

p(g, T) = e- lHT A T ( P0 )e iHT = e~^ T A T ( P o)e lh9T , (2.1) 

where At is the superoperator that describes the cumulative effect of decoherence in the channel. 
The final state can always be written in the form (|2.ip . by going to an interaction picture relative 
to the Hamiltonian (jl.2p . but in doing so, the decoherence superoperator, At, generally becomes 
dependent on g. In our analysis, however, we assume all decoherence processes to be invariant under 
rotations about the z axis (we also assume that the decoherence is independent and identical from 
one qubit to the next), which implies that At is independent of g and also means that At commutes 
with e~ %hgT , i.e., p(g,T) = AT{e~ thgT poe th9T )- The final state contains the information about the 
value of g. The accuracy with which we can distinguish the state p(g, T) from neighboring states 
on the one-parameter path parametrized by g controls the accuracy in the estimate of g. 

From the results of measurements on a set of v probes, we obtain an estimate g e st for the value 
of g. We can quantify the statistical deviation of the estimate from the true value of g by the 
units-corrected deviation from the actual value, 

Htraw-'Jf ■ (2 ' 2) 

introduced in . A lower bound on 5g is given by the generalized uncertainty relation 

89 > 1 , (2.3) 
y/v {ds/dg) 



where ds is the "statistical distance" between neighboring quantum states along the trajectory 
parametrized by g. The statistical distance is given in terms of the change, 

dp = ^^dg ee p'dg , (2.4) 
dg 



in p = p(g, T) due to a small change dg in the value of g: 

~{p'C p (p')) . (2.5) 



ds\ 
dg) 



2 

= trl 



In the basis {\a}} that diagonalizes p(g,T) = J2 a Pa\ a )( a \i the superoperator C p takes the form 

C p (0)= ]T -^—O aP \a){f3\. (2.6) 
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The operator C p (p') is called the symmetric logarithmic derivative [2j because 

pi= l -(pC p (p')+C p {pi)p) . (2.7) 

The quantity (ds/dg) 2 = ti[p' C p (p')] is often called the quantum Fisher information 0]. 

The generalized uncertainty relations are derived using a quantum version of the Cramer-Rao 
bound p|, 0, 0]- Generally, this bound can be achieved only in the case of an optimal measurement 
on each probe and, even then, only asymptotically for a large number v of probes, as emphasized 



by Braunstein [271 ] . In our analysis, we explicitly exhibit an optimal measurement, and we let f m j n 
denote the number of probes required to approach the bound within some fixed fractional error. 

Consider now the the continuous path in the space of states of the probe parametrized by g. 
Nearby points on this path are related by the derivative 

p' = -iT[h,p] = -iT[h,p] =iT^2(p a - pp)h a p \a){/3\ , (2.8) 

where here and throughout a hat denotes the difference between a quantity and its mean value, 
i.e., h = h — (h). Plugged into Eq. (|2.6p . this gives 

C p (p') = 2iT^2^^-h aP \a)(P\ (2.9) 



and 

" ^ 4T 2 A 2 < AT 2 {Ah) 2 , (2.10) 



dgj 



where 

\2 



2 t5 p " + 1 'i i 



and 



(Ah) 2 = (h 2 ) = l Y,^+Pp)\Kp\ 2 (2.12) 

a,f3 

is the variance of h with respect to p(g,T). Notice that in Eqs. (|2.9ft and (|2.1ip . we can remove 
the hat from h without changing anything, whereas in the variance (|2.12p . we cannot do so. From 
Eqs. &3$) and (l2~T0"l) . we obtain 

59 2 V^fA 2 T^TAk ■ (2 ' 13) 

These are the generalized uncertainty relations in the forms we will use. Notice that when p(g, T) is 
a pure state, we have A = Ah, and thus equality holds in Eq. (|2.10p and in the second inequality in 
Eq. (|2.13p . Notice also that in the absence of decoherence (i.e., when At is the unit superoperator), 
a pure initial state stays pure, and Ah is independent of T. 
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The second (weaker) inequality in Eq. (|2. 13|) is generally easier to work with than the first 
(stronger) inequality, because computing the uncertainty Ah is usually easier than computing the 
corresponding term in the first inequality. In the case of an initial pure state in the absence of 
decoherence, the two inequalities are equivalent. In this case, to minimize the uncertainty in g, 
we should initialize the probe in a state in which Ah is maximal. We show in Sec. IIVI that when 
there is decoherence in the probe qubits, which changes an initial pure state to a mixed state, the 
weaker inequality is not very useful, and we are forced to work with the stronger inequality. 

For nonGaussian statistics, the quantum Cramer-Rao bound is saturated only in the limit 
v — > oo, i.e., when the measurement process involving an n-qubit probe is repeated many times. 
We let 

^min denote the minimum number of iterations that are required for the measurement 
accuracy to approach the quantum Cramer-Rao bound within some fixed fractional error. For the 
qubit protocols we analyze, v m \ n is essentially independent of protocol, as we discuss further below 
when we consider measurements that achieve the bound f)2. 13|) : a typical value might be, say, 50. 
The need to do at least f m i n iterations places a constraint, v > f m in- Together with the constraint 
that each probe must contain at least one qubit, this gives us the following constraints on n: 

1 < n = N/v < iV/f m in • (2.14) 

For these constraints to be consistent, it must be true that 

R(t — T) = N > i/ min . (2.15) 

As we discussed in the Introduction, we assume in our analysis that the parameter must be 
determined in a fixed time r = vn/R + T and that the qubit supply rate R is a fixed resource. For 
a particular kind of decoherence, we vary the interaction time T and the number of qubits in each 
probe, n = R(r — T)/v, to achieve the best accuracy in determining g. 

A different sort of resource that is dependent on the way the probe is initialized is the magnitude 
of A. If we assume that the energy spread for each of the qubits is fixed, then a way of getting a 
large value for A is to initialize the n qubits in each probe in an appropriate entangled state. In a 
sense, A is itself a measure of the entanglement or quantum coherence available for improving the 
ability to determine g. 

III. MEASUREMENT ACCURACY IN THE ABSENCE OF DECOHERENCE 

In this section we review the limits on the accuracy of estimating g in the case where there is no 
decoherence. From the generalized uncertainty relation, we see that the initial state of the probe 
has a direct bearing on the optimal accuracy. We look at two very different initial pure states 
of the probe. In the first case the n probe qubits are initialized in a product state, and in the 
second case they are in a collective entangled state. For both cases we compute the limit on the 
precision with which g can be estimated. Since there is no decoherence, the probe state remains 
pure, and the two inequalities in Eq. f)2 . 13|) are equivalent, because A = Ah. Thus in this section 
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we only need to consider Ah. This section also serves to establish our notation and to summarize 
the results in Q]. 

A. Initial pure product state 

If the probe is initialized in a pure product state, p p , of the n qubits, we have 



dsp _ 

dg \\ j± x \dg J 

as shown in Appendix [A] (see also [6]). Here and in the following, the subscript p stands for 
"product state." The line element dsp is in the space of n-qubit density operators, while dsj are 
line elements in the space of states of each the n probe qubits. From Eq. (|2.13p we see that the 
best choice of initial probe state is one that maximizes 



A = Ah 



\ 



E( A ^) 2 > ( 3 - 2 ) 

3=1 



where (Ahj) 2 is the variance of hj for the jth qubit. Thus we have to maximize Ahj for each of 
the n qubits, and we do so by initializing each of the n qubits in a pure state lying in the equatorial 
plane of the Bloch sphere of states for each qubit. Here we choose initial state 

l^> = -^(|0i) + |li» or Pj = l(l j + (Tx;j ) . (3.3) 

The vectors \0j) and denote the eigenstates of a z -j for the jth qubit. The initial state of the 
probe is 

n 1 n 

3=1 3=1 

The effect of the coupling to the parameter is to rotate the Bloch vector of each of the qubits 
around the a z axis. At time T, after passage through the channel, the state of each qubit has 
rotated through an angle gT, giving a probe state 

n 1 n 

P P (g, T) = Pj (g, T) = — (g) (lj + a x , 3 cos gT + a y .j singT) . (3.5) 

3=1 3=1 

We use the evolved state (|3.5[) in our discussion of achieving the optimal measurement accuracy 
in Sec. IIII CI For the present, however, since the variance of hj is unchanged by the evolution, we 
can evaluate it using the initial state (|3.3p . This gives a variance (Ahj) 2 = 1/4 for each qubit, and 
thus 

(Ah) 2 = n(Ahj) 2 = j . (3.6) 
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The generalized-uncertainty bound on the estimate of g becomes 

6g - dgp= T^ = WW = T^R(T-T) > (3 - 7) 
where N is the total number of qubits used in the measurement scheme, and 5g p is the bound for 
pure product inputs. 

Since the bound depends only on N, and not separately on n and v, we can always choose n = 1 
without affecting the optimal measurement accuracy. This, of course, is the statement that for 
product-state inputs, the measurement accuracy is indifferent to whether we regard the qubits as 
gathered together into multi-qubit probes. To find the optimal bound, all that is left is to adjust 
the interaction time T to minimize 5g p . Doing so gives 

T = \r (3.8) 

and thus 

N = v = ^Rt . (3.9) 

That there be enough probes to satisfy v > u m - m requires that \Rt > u m i n . When ^Rt < v m i n , 
the measurement bound is optimized by the choices n = 1 and v = t> m in- This gives an interaction 
time T = t — v m i n /R that decreases with r until Rt = v m i n , at which point it is impossible to 
obtain and use i/ mm probes within the overall duration r. This interaction time occurs in every 
situation we consider, when the measurement protocol is starved of qubits, so we abbreviate it as 

T s = T-v min /R. (3.10) 

The optimal bound on measurement accuracy thus takes the form 

1 



1< Rt < 3v, 



nun > 



S 9P = < 3 ^3 /2 (3.11) 
3z^ min < Rt. 



tVRt 

In our resource-based analysis, in which the overall measurement time r and the rate R at which 



qubits can be supplied are the resources, the 1/V Rt 3 scaling is the signature of the standard 
quantum limit. The behavior of the bound for 1 < Rt < 3i/ m j n is included for completeness in our 
subsequent analysis, but is not so important since it expresses what happens when the measurement 
protocol is starved of qubits. 

B. Initial pure entangled state 

If the probe can be initialized in an entangled state, we can obtain bigger values of Ah. The 
maximum value is obtained by superposing two n-qubit eigenstates of h corresponding to the lowest 
and highest eigenvalues. Thus we initialize the probe in the "cat" state 0] 

|*c) = -^(|00...0) + , (3.12) 
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denoted by the subscript c. The initial density operator p c = |^ r c )(^ r c | can be written in the form 

^ / n n n n \ 

^=2^+T ® + ^ + ® " ^aO + ® ^ + + ® " **vu) • ( 3 - 13 ) 
\j=i i=i i=i j=i J 

After passage through the quantum channel, the state of the probe becomes 
Pc(g,T) = 

* ' \i=i i=i 

n n \ 

-i^) . (3.14) 

3=1 J 

We use this form in our discussion of achievability in Sec. IIIICI Since Ah does not change under 
the quantum evolution, we can evaluate it using the initial cat state, which gives 

(Ah) 2 = ^ . (3.15) 
The generalized-uncertainty bound becomes 

Sg > Sg c = 1 = !— = . (3.16) 

Tn^v Ty/Nn TR(r - T) 

By letting the probe qubits be in a collective entangled state, the accuracy in our estimate is 

enhanced by a factor of 1/ y/n. 

For the cat-state input with no decoherence, it is optimal to make v as small as possible, i.e., 

v = fmm. This puts as many qubits as possible into each probe consistent with the constraint ([2.140 . 

i.e., n = N/v m i n , which is clearly optimal in the absence of any decoherence to degrade the 

entanglement. To find the optimal bound on the measurement accuracy, we adjust the interaction 

time T to minimize Sg c , giving 

T = \r (3.17) 



and 



N = nv min = X -Rt . (3.18) 



In order that v > v m i n , we require \Rt > iy m - m . When ^Rt < f TO i n , we choose n = 1 and u = u n 
which gives T = T S and an optimal bound that is the same as for product-state inputs. 
The optimal bound on measurement accuracy thus becomes 

' 1 < Rt < 2u min , 



6g c = { ^)/^ (3 19) 

Rt2 , Rt > 2v min . 

The 1/Rt 2 scaling of the cat-state bound is the signature of the so-called Heisenberg limit; it 
is to be contrasted with the corresponding l/V Rt 3 scaling available from product states. The 
enhancement available from entanglement is roughly a factor of \/^min/R T - Just as for product 
states, the behavior of the bound for 1 < Rt < 2v m \ n is not so important, as it expresses what 
happens when the measurement protocol is starved of qubits. 
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C. Achieving the optimal measurement accuracy 

Estimating g involves making measurements on the probe qubits. A strategy that gives an 
optimal estimate of g is to measure a x on all the probe qubits. This provides a means of estimating 
gT, from which g can be calculated provided we know the interaction time T accurately. 

In the case of an initial pure product state, we can specialize to having just one qubit in each 
probe (n = 1), prepared in the state f)3. 3|) . We measure a x on each qubit. The results are averaged 
over v trials to get an accurate estimate for g [9]. The expectation value and variance of o x with 
respect to the evolved state in Eq. (|3.5p are 

(a x ) = cos gT , (Aa x ) 2 = 1 - cos 2 gT = sm 2 gT . (3.20) 

The average of the results over v trials, which we denote a x , has the same expectation value, but 
its variance decreases by a statistical factor of i.e., Aa x = \ sin gT\ / ^/i' . 

We estimate g as g es t = T -1 arccos a x . When the uncertainty in a x is small enough, we can 
approximate (g cs t) = T arccos (a x ) = g and 

5g = Ag cst = A ^ = . (3.21) 

\a{<T x )/dg\ TyV 

The approximation here is that the datum a x must be likely to lie close enough to the expected 
value (a x ) that a linear approximation to the arccos function at the operating point is valid. This 
requires that v be large enough that 1 3> Aa x = Aa x /^Jv ~ \j*Jv. That v must be large is the 
expression, in the context of this particular measurement, of the general fact that the quantum 
Cramer-Rao can only be achieved asymptotically; it leads to our requirement that v > z^ m i n S> 1. 

When the probe is initialized in the cat state, an optimal measurement strategy is to measure 
o x on all n qubits simultaneously and to multiply all the results together Qj. Formally, this 
corresponds to measuring 

■n 

S x = (g)(T x;J (3.22) 

3=1 

The expectation value and variance of T, x with respect to the evolved state (|3.14|) are 

(E x ) = cosngT , (AT> X ) 2 = 1 — cos 2 ngT = sin 2 ngT . (3.23) 

The average of the results over v probes has the same expectation value, but its variance decreases 
by a statistical factor of 1/u. 

We estimate g in the same way as above for product inputs. The only difference is the additional 
factor of n in the rotation angle due to the coherent rotation of the entangled qubits in each probe. 
The resulting uncertainty in our estimate of g is 

5g = — L= , (3.24) 

1 Jly/V 
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thus saturating the bound (|3. 16|) . Notice that the condition for making a linear approximation to 
the arccos function is the same as for product inputs, i.e., 1 ^> AT, X = AT, x /y/u ~ showing 
that we can take z^ m i n to have the same value for product and cat-state protocols. 

There are technical questions associated with how one resolves the fringes in Eqs. (|3.20|) and 
(|3.23p in order to zero in on the actual value of g. These questions are well understood, however, 
and are irrelevant to our goal of understanding the effects of decoherence, so we do not consider 
them further. 



IV. MEASUREMENT ACCURACY IN THE PRESENCE OF DECOHERENCE 



The previous section reviewed, within the context of our resource-based analysis, the measure- 
ment accuracies that can be obtained in the absence of decoherence. In this section we introduce 
decoherence to see how it affects the accuracy of parameter estimation. We consider a general 
model for decoherence of the probe qubits, subject to the restrictions that the decoherence (i) is 
independent and identical from one probe qubit to the next, (ii) is continuously differentiable and 
time stationary, and (iii) commutes with rotations about the a z axis. Since the interaction Hamil- 
tonian that connects the probe qubits to the parameter generates rotations about the a z axis of 
each of the qubits, the effect of the third restriction is to separate cleanly the effect of the parameter 
from the effects of decoherence. 

Decoherence can be described in terms of trace-preserving quantum operations (completely 



positive maps) on density operators [28|, |29|, |3C], [3l|, [32J, [33J . A quantum operation on single-qubit 
states is completely specified by the transformations of the operator basis set consisting of 1, a x , a y , 
and a z . A general time-dependent trace-preserving map At on one-qubit states, which commutes 
with rotations about the a z axis, has the form 

A t (l) = t + f(t)a z , 
A t (a z ) = g(t)a z , 
At (cr x ±i(T y ) = h±(t)(a x ± iay) , (4.1) 

where f(t), g(t), and /i+(t) = h*_(t) are arbitrary functions of time t. The requirement that 
the evolution described by At be continuously differentiable and time stationary implies that the 
derivatives of the quantities on the left of Eq. (|4.ip be linear combinations with constant coefficients 
of these same quantities. Thus we have 



dA t (l) 
dt 

dAt(a z ) 
dt 

dAt(cr x ± icr y ) 
dt 

where fi, 71, 72, and uj are real constants. 



HliA t {<Jz) , 
-liAt{a z ) , 

-(72 ± iuo)A t {(J x ± i(T y ) , (4.2) 
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The solution of Eqs. (|4.2p . with .Ao = gives the most general single-qubit decoherence model 
that satisfies restrictions (ii) and (iii) above: 

A t (l) = 1 + M(l " e~ 7l > z , 

A t {a x ± iff,,) = e^e^\a x ± iff„) . (4.3) 

The dissipation in this model is that of the standard qubit decoherence model, involving a longitu- 
dinal decay time T\ = I/71 and a transverse dephasing time T2 = 1/72- In the limit i — * 00, every 
single-qubit state decays to the state i(l + /xcj^), which means that we must have — 1 < /x < 1. 
Complete positivity requires that T2 < 2T\. 

In addition to the dissipation, there is a coherent rotation about the a z axis by an angle cot. 
If a decoherence process does introduce such a coherent rotation, it cannot be distinguished from 
the rotation produced by the parameter g; any procedure for estimating g would actually estimate 
g + to. Throughout the following, we assume that the decoherence model does not include any 
coherent rotation, but for convenience, we incorporate the rotation due to g into At by assuming 
that ijj = g and omitting the further coherent rotation in Eq. (|2.ip . The mapping of the Bloch 
sphere induced by At (with u = g) is illustrated in Fig. [TJ 




FIG. 1: (color online) The transformation of the Bloch sphere under the map At of Eq. (|4.3[) at 72^ = 0, 
7r/4 and 7r/2, with /1 = 3/4 and 71 = 372/4. The dotted lines show the unit sphere. The coherent rotation 
due to lj = g cannot be seen in these diagrams. 



A. Initial pure product state 



We are now prepared to investigate how decoherence affects the theoretical minimum for 5g 
when the probe is initialized in a product state. At time t after entering the quantum channel, the 
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state of a probe is given by applying the map (|4.3p to the initial state of each qubit in the pure 
product state (|3.4p : 

1 n 

Pp(9, t) = ^ (8) (% + M 1 - e" 71 *)cr 2;i + e" 72t {a x .j cos gt + sin^)J . (4.4) 
i=i 

We first look at the weaker inequality in Eq. (|2.13p . At the time t = T when the probe leaves the 
quantum channel, we have (hj)x = ^ti[p p (g, T)a z -j] = p(l — e~ 7lT )/2 and (^|)t = 1/4, giving a 
variance 

(A/ l ) 2 = ^) T = ^(l- M 2 (l-e-^) 2 ) . (4.5) 
The resulting weaker uncertainty-principle bound from Eq. (|2.13l) is thus 

S 9 > SgtoM = 1 1 . (4.6) 

Ty/N ^ _ ^2 (1 _ e -7iT)2 

Decoherence in the transverse {^x;j-^y;j) plane does not appear explicitly in the bound 6gp (71). 
This is because the weaker inequality in Eq. (|2.13|) is determined by the variance of h, which 
depends only on the decoherence along the longitudinal a z -j direction. 

Our conclusion is that we should not rely on the weaker inequality in Eq. (|2.13p to provide a 
good bound on the maximum achievable measurement accuracy when there is decoherence. For 
instance, in the case where there is only transverse decoherence, i.e., 71 = 0, the bound 5g p w (71) 
remains constant at l/Ty/N, even though the transverse decoherence ultimately leaves the probe 
qubits in a state along the a z axis where the rotation produced by the parameter has no effect. To 
see the dependence of the measurement accuracy on the transverse decoherence, we have to use 
the stronger inequality in Eq. (|2. 13|) . 

Turning to that stronger inequality, we need to evaluate A as in Eq. (|2,lip . and for that purpose, 
we first write the state of each probe qubit after passage through the channel in diagonal form, 



where 



and 



Pj (g,T) = l(l+ sjd* + dA l+X+l + 1(1- y ^ + 4 j l-X-l < 1-7) 



di = fjt{l ~ e" 7lT ) , d 2 = e~ 72T (U 



|+) = cos(0/2)|O) + e i9T sm{9/2)\\) , 

|-) = sin(0/2)|O) -e^ T cos(0/2)|l) , (4.9) 



are the eigenstates of pj(g,T), with 



sin 9 = 



d 2 

\M 2 + di 



(4.10) 
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To evaluate Aj for the jth qubit, we need the off-diagonal matrix element of hj = \& z -j in this 
eigenbasis: 



(hj)+- = (hj)!+ = 2<+|er,y| 



ifrini 



Plugging Eqs. (|4"7T|) and (|4"TTj) into Eq. (|2"TT|) . we find that 



A'- 



nA 2 = -d 2 2 

3 4 z 



™ -2 72 T 

4 



from which follows a stronger uncertainty-principle bound, 



(4.11) 



(4.12) 



(4.13) 



TVN T^R(t - T) ' 

for v uses of the quantum probe. This is a much more reasonable bound on measurement accuracy, 
since it depends explicitly on the transverse decoherence that we expect to make a difference in the 
measurement; moreover, the dependence simply degrades the measurement accuracy exponentially 
with the number of T 2 times for which each probe is in the quantum channel. 

We can show that directly that the bound (|4.6p is weaker than that of Eq. (|4,13p through the 
following chain of inequalities: 



1 - /x 2 (l - e~^ T ) 2 > ^1 - (1 - e~^ T ) 2 = ^e~^ T {2 - e~^ T ) > e" 7lT/2 > e 



-72 T 



(4.14) 



the last of which requires the complete-positivity condition, 72 > 71/2. 

Since the probe state is separable at all times, we can choose n = 1 without affecting the optimal 
measurement. What remains is to choose the interaction time T, within the range < T < r, 
so as to minimize the bound 6g^\j2)- There is a single minimum at T = T p , determined by the 
equation (72 T p ) 2 — (3/2 + 72t)72T p + 72T = to occur at 



3/2 + 72 T - ^(3/2 + 12 t) 2 - 4 72 r 



(4.15) 



We cannot use this interaction time when it becomes so short that the measurement protocol is 
starved of qubits, i.e., when R{j — T p ) < v m in- in this situation, we choose n = 1 and v = f m i n , 
which gives the interaction time T s of Eq. (|3.1U|) . 

Plugged into Eq. (|4.13j) . these interaction times give the optimal value of the bound Sgp (72)) 
which can be written in the dimensionless form 



72 



=72 T B 



72 

e 72T p 



^ , I 'm 1 ^min \ 

< 72-r < 72 [Tp + —— J . 



l2T pX h 2 (r - T p ) 



R - ' "\-- ' R 
72T > -f 2 [T p + ——j. 



(4.16) 



We plot the dimensionless optimal interaction time, 72 Tp, and the resulting dimensionless optimal 
bound (|4.16p as functions of 72T in Fig. [2 

There are two important limits. When the transverse decoherence has little effect during the 
overall time r, i.e., 72T < 1, we find T p = 2r/3 and an optimal bound that reduces to the optimal 
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FIG. 2: Dimcnsionless optimal interaction time, 72 T P of Eq. (|4. 1 5[> . and dimensionlcss optimal bound (|4 . 10[) . 
plotted as functions of dimensionless interaction time 72T. In these plots we assume that R is big enough 
that we do not encounter the situation where the protocol is starved of qubits, since this situation is of little 
interest. 

measurement accuracy (|3,lip in the absence of decoherence. In contrast, for large transverse 
decoherence, i.e., 72T 3> 1, the optimal interaction time is T p = T2, and the optimal bound becomes 
$9p (72) = e/T2\fRr or, in terms of the dimensionless optimal bound, \J R/^bg^ H2 = e/^/72T. 
In this case, it is optimal to have each qubit scoot through the quantum channel in a dephasing time, 
before the dephasing can destroy the effect of the parameter-induced rotation; roughly speaking, 
each qubit determines g with accuracy e/Ti, which is improved by the statistical factor 1/y/Rr 
corresponding to the number of qubits used in time r. 



B. Initial pure entangled state 



We now look at the case in which the probe is initialized in an entangled cat state. The density 
operator of the probe, after a time t in the channel, is obtained by applying the map (|4.3|) to the 
initial cat state f)3. 13|) : 

1 / n n 

Pc(9, t) = <g) (l, + [e-»* + Ml - e~ 71 ')K;i) + ® (*j " [^ 7l * - ^ ~ e~^)]a z , 

\i =1 J =1 



3=1 3=1 



(4.17) 



To evaluate the weaker inequality in Eq. (|2.13p . we need to evaluate the variance of h. Using 
Eqs. (jBip . (|B2p and (|B3[) from Appendix [Bl we find that at the time t = T when the probe exits 
the quantum channel, (h)x = nd\/2 and 

(h*) T = \{l + {n-l){e-** T + ti)). (4.18) 
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From these we obtain the variance of h at the time T when the probe exits the quantum channel: 

(Ah) 2 = (h 2 ) T = j (l + (n - l)e~ 2 ^ - dfj . (4.19) 
The weaker bound in Eq. (12. 13ft on the accuracy of estimating g thus becomes 

Sg > SgtHli) = T \ N 1 1 • (4.20) 

T ^ N Jl + (n- l)e- 2 ^* -ii 2 [l- e-~nt f 



We do not expect this bound on 8g c {l) to be particularly useful because, as for the case of an 
initial product state, the decoherence in the transverse directions for the qubits does not come into 
the bound at all. 

We now look at the bound on measurement accuracy placed by the stronger inequality in 
Eq. (|2.13p . The density operator p c (g,t) is diagonal in the tensor-product basis formed by the 
eigenvectors of cr z .j, except in the two-dimensional subspace spanned by the vectors 1 00 . . . 0) = |0) 
and 1 1 1 ... 1) = |1). We denote this subspace of the n-qubit Hilbert space TL by /C. The operator 
h is diagonal in the tensor-product basis formed by eigenvectors of u z -a. From Eq. (|2.1ip . we see 
that there is no contribution to A 2 from the subspace in which p c (g,t) and h are simultaneously 
diagonal. Thus, for computing A 2 , we can work with the operators p c (g, t) and h that are obtained 
by projecting p c (g,t) and h down to the subspace fC, i.e., 



A 

and 
with 



.(g,T) = i (d+|0)(0| +(f 2 n e- m£,T |0)(l| + d?e ingT \l}(0\) (4.21) 



S = ^(|0>(0|-|1>(1|) , (4.22) 



d±s (i±£^)" + (i^^)\ ( , 23 ) 

The next step is to write p c (g,T) in diagonal form, 

p c ( g ,T)=p + \+)(+\+p_\-)(-\ , (4.24) 

where 

1 



P± = I [d+ + d_ ± y/(d+ - d_) 2 + 4d^j (4.25) 
are the eigenvalues of p c and 

|+) = cos(0/2)|O) + e in9T sin(#/2)|l) , 

|-) = sin(0/2)|O) -e insT cos(0/2)|l) , (4.26) 

are the eigenstates, with 

sin = m i . (4.27) 
1 {d + - d_) 2 + Ad 2n 



16 



To evaluate A, we need the off-diagonal matrix element of h in this eigenbasis: 

— — — n 

h+- = h*_ + = (+\h\-) = ^mnd. (4.28) 



The resulting value of A is 



{P + - P-?\h + -\ 2 = ^e~^ T , (4.29) 



from which follows the stronger uncertainty-principle bound for a cat-state input, 

e n 72 T n l2 T r~ 12 RT{t-T)/u 

Aside from being stronger than the bound (|4.20|) . this is a more sensible bound, since it depends 
explicitly on the transverse decoherence. When 72 = 0, this bound simplifies to the cat-state bound 
in the absence of decoherence, Eq. (|3,16p . Moreover, by comparing with the bound for a product 
input, Eq. (|4.13p . one sees that this bound retains the 1/y/n advantage purchased by using an 
entangled input, but at the price of a decoherence rate that is n times faster. 

We have assumed that both n and T are controllable parameters in the estimation scheme we 
are considering, with v determined by Eq. (jl.4p . To minimize the bound (|4.20p . we use the second 
form, from which v has been eliminated. The values for n and T that minimize can then be 
found by solving simultaneously the two equations, 

dSg { c s) e n ^ T 

= = (2ra 72 T - 1) , 

dn 2nTy/nR(r - T) V ' ' 

= = ; [3T + 2n 72 r(r - T) - 2r] , (4.31) 

which give n = I/72T and T = r/2. The determinant and trace of the Hessian of 5gi s \ with respect 
to n and T, evaluated at this point, are both positive, showing that it is indeed a minimum. The 
minimum value of the bound is 

6gl'){v) = -^=. (4.32) 

This minimum cannot always be attained, however, because we have the additional constraints 
of Eq. (|2.14p . i.e., 1 < n < N/v m - m = R(r — T)/u min = RT/2v m - m , which do not always allow us to 
choose n equal to the optimal value 1/727"- If 72 T does not lie between 2v ra \ n /Rr and 1, we have to 
choose a value for n that lies on the boundary of allowed values. There are two cases to consider. 
If the decoherence rate is high, i.e., 72T > 1, we choose n = 1, thus using probes consisting 
of individual qubits to estimate g, in which case the analysis reduces to that of the preceding 
subsection. Notice that 72T = 1 gives j2(T p + v m m/R) = 1/2 + ^ m in72/-R- Thus if 2v m \ a r y2/R < 1, 
the second case in Eq. (|4. 16[) applies whenever 72T > 1. If, however, 2v m ;ml2/R > T the protocol 
begins to be starved of qubits for some 72T > 1, and there is no situation where cat states offer 
any advantage. Throughout the following, therefore, we assume that 2v m \ n ^f2/R < 1- 
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TABLE I: Minimum value of bound on estimating g using cat-state probes with available resources deployed 
optimally. The table assumes that 2i/ m i n 72/i? < 1. 

Range of 72 t T n v N Sgi s (72) 

< 72T < 5 — / - of Eq. i l.i.lUl i 1 ;/„,,,. - 



t-» : a • , ~\ a wi ■*— m • ' » V l " 111111 ** 111111 m > 



2f mi n72 ^ . /2f raill 72 _ /n d_ /<-,.. .. d_ /o 4 V^min p72 flr 2 /4^ 



2^ mi n72 , „ , _ /n , _ .. D _ 2/0 D _ /0 2V2e 



< 7 2 r < 1 r/2 1/72 r 7 2 i?r 2 /2 i?r/2 



11 Ty/R/j 2 

7 2 r > 1 T p of Eq. (@35]) 1 i?(r - Tp) i?(r - T p ) 



e 1 



Tpy^r - T p ) 



If the decoherence is small, i.e., 

72 r < 2^ min /i?r ^ r 2 < 2^ min / 72 i? , (4.33) 

we use the largest cat state that can be constructed from the available resources, thus choosing 
v = v m - m . Using the last form in Eq. (|4.30p . with v = i^ m i n , we find that bg£^ has extrema for 
T(t — T) = and T = r/2. We discard the first possibility because it is inconsistent 

with the constraint (|4.33j) . i.e., ^min/72-R = T{r — T)< r 2 /4 < u m - m /2^/2R- Moreover, the second 
derivative of Sg^ with respect to T, evaluated at T = r/2, is strictly positive when Eq. (|4.33p is 
satisfied, showing that T = r/2 gives a minimum. The optimal interaction time is again T = r/2 
(n = i?r/2 

^min)) and the minimum value of the bound becomes 

§g (s) ( 72 ) = M^ e 7 2 ^ 2 /4*W (4 34) 

Notice that Eq. (|4.34|) reduces to the second case in Eq. (|3.19p when there is no decoherence, i.e., 
when 72 = 0. 

The first case in Eq. (|3.19p reminds us that one further case occurs at very short times, when 
the protocol is starved of qubits. In particular, when Rr/2v m \ n < 1, we must choose n = 1 and 
v = i^min, leading to the familiar interaction time T s of Eq. (|3.10p . 

We can now piece together the various regions that govern the optimal bound on the estimate of 
g using cat-state probes with the available resources deployed in the optimal fashion. The results 
are summarized in Table HI The top row is the case where the protocol is starved of qubits; the 
second row is the case of low decoherence, for which the probes are prepared in cat states containing 
as many qubits as allowed by the need to have at least f m i n probes; the bottom row is the case 
of high decoherence, for which the probes are individual qubits; and the middle row describes the 
transition from high decoherence to low decoherence. The two regions where cat states play a role 
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VRSg^ p (y 2 )/y 3 2 




1.x10 -6 0.0001 0.01 1 

FIG. 3: (color online) The four thick lines labeled (a), (b), (c), and (d) show the dimensionless optimal 
bound (|4.35p for cat-state inputs, plotted as a function of dimensionless interaction time 72T, for v m i n = 50 
and y/R/2v min j 2 = 1, 10, 100, and 1 000, i.e., y/R/-y 2 = 10, 100, 1 000, and 10 000, respectively. Note that 
both axes use a logarithmic scale. The use of cat states provides no advantage for 72T > 1. The two regions 
where cat states provide an advantage (second and third rows of Table |TJ are absent for ^/i?/2i/ m ; n 72 = 1, 
but become apparent for the other three values of y/ -R/2i/ m i„72. In terms of the dimensionless bound, 
the transition region between high and low decoherence (third row in Table [TJ) has a form independent 
of -R/72, but extends to smaller values of 72T as -R/72 increases. For small enough 72T, the protocol is 
starved of qubits, and cat states again provide no advantage over product states. The thin dotted lines show 
the dimensionless product-state bound (|4.16p for the same four values of yj -R/2^ m in72. The product-state 
bound agrees with the cat-state bound in the high-decoherence region (72 t > 1) and with the corresponding 
cat-state bound in the region where the protocol is starved of qubits (f m in72/-R < 72^ < 2f m i n 72/-R); in 
between, where cat states provide an advantage, the product-state bound is independent of -R/72. 



exist when 2^ m j n 72/i? < 1. The dimensionless bound introduced in Eq. (|4.16p is given by 



72 72 



=72 r s 



72 72^-vAw ' 



^min72 , 2f min 72 
< 72T < 



R ( 72 r) 2 



min e (R/72 ) (72 r ) 2 /4i/ min 



R 

2^min72 

R 



R 



< 72T < 



2f m in72 



R 



12T 



=72^ 



2^min72 

R 

7 2 t > 1. 



(4.35) 



< 72T < 1, 



_ T^T p ^72(t - T p ) 

This dimensionless optimal bound is plotted in Fig. [3] for the choice v n 



50 and for several 



values of \JRjyi- 

The conclusion to be reached from the results summarized in Table [J and Fig. [3] is that cat-state 
entanglement is only useful for improving the estimate of g when one wants to estimate g on a 
time scale that is shorter than the time scale over which decoherence acts. 
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C. Achieving the optimal measurement accuracy 

We now examine the effectiveness of the measurement strategies described in Sec. 1111 CI in the 
presence of decoherence in the probe qubits. When the input state is a product state, we can again 
specialize to the case where there is only one qubit in each probe (n = 1) and v = N. We measure 
o x on each qubit and average over the results of v trials to obtain the estimate of g. After an 
interaction time t = T, the expectation value and variance of a x for a single qubit in the evolved 
state (IP1) are 

(a x ) = e" 72T cos gT , (A^) 2 = 1 - e" 272T cos 2 'gT . (4.36) 

Averaging over the results of v trials yields the quantity a x , which has the same expectation 

value as a x , but has variance reduced to (Aa x ) 2 = (Aa x ) 2 /v, and we estimate g as g = g es t = 

arccos(e 72T o' a ;). After many trials, the variance of a x becomes small enough that we can 

approximate (gest) = T~~ l arccos(e 72T (<7 x )) = g and 

A Aa a . e 7 2 T ,/i _ e -2 72 r cos 2 gT 

5g = Ag est = = —= — — — . 4.37 

\d(a x )/dg\ T^/v |smsT| 

If we use this straightforward method of estimating g by averaging measurements of a x on all qubits, 

achieving the bound (|4.13p for determining g, with the interaction time adjusted to the optimal 

value T = T p , requires |singT p | = 1. Even though g is not known, this can be accomplished 

by using feedback onto the rotation of the qubits to find this sweet spot on the fringe pattern. 

An alternative to feeding back onto the rotation of the qubits is feedback to rotate the quantity 

measured in the equatorial plane of the Bloch sphere until the desired operating point is achieved. 

When the probe is initialized in a cat state, we have seen that the optimal strategy, in the 

absence of decoherence, is to measure a x on all n qubits simultaneously. The same measurement 

works when there is decoherence. Using Eqs. (|3.22p and (|4.17p . we find that after an interaction 

time t = T, 

(Sju) = e" n72T cos ngT , (AS a .) 2 = 1 - e" 2n72r cos 2 n 5 T . (4.38) 

We can convert the results for product-state inputs to this case by the substitution T — > nT, so 
the average over v probes leads to 



e n72T y/1 - e ~ 2n ^ T cos 2 ngT lA . 

bg = - — -= — — . 4.39 

1 riyjv | suing 1 \ 

We can saturate the bound (|4.30p on estimating g, for the appropriate interaction time T, by using 
feedback to operate at a point where | sinngT| = 1. 

One last point concerns the question of making a linear approximation to the arccos function, 
which allows us to relate the mean and variance of <7 es t directly to the mean and variance of S x . 
This approximation requires that v be large enough that 1 3> e n ^ 2T AS x = e™ 72T AS x / \Jv. The 
results summarized in Table |T] show that for the optimal choices of n and T, it is always true that 
7272T ~ 1, showing that the requirement is that v be large in a way that is independent of the 
details of the protocol. 
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V. DISCUSSION 

In this paper we have studied quantum limits on determining a frequency g that controls the 
rate at which qubits rotate about the z axis of the Bloch sphere. The question of determining 
g is the same as the problem of distinguishing qubit states that differ by having been subjected 
to different rotations. The quantum limits on determining g are well known [9]: if n qubits are 
prepared in product states, the uncertainty in determining g scales as 1/y/n, the standard quantum 
limit or shot-noise limit, whereas if the same n qubits are prepared in an entangled cat state, the 
uncertainty scales as 1/n, which is called the Heisenberg limit. Our purpose in this paper has 



been to investigate, following [10j], how these scalings change when the qubits are subjected to 
independent, but identical decoherence processes during the measurement. 

The decoherence model we consider is the most general continuous-time process that is invariant 
under rotations about the z axis. This results in a standard qubit decoherence model in which 
qubits decay with a time constant T± and lose phase coherence with a time constant Ti- It is the 
phase-coherence time Ti that is important for efforts to estimate g; cat-state entanglement is only 
useful for times of order or smaller than T2 . 

To make the analysis meaningful, we introduce as resources the rate R at which qubits are 
supplied and the overall time r that one has available for estimating g. If one does not introduce 
these resources, one can achieve any desired accuracy in estimating g by assuming that one can 
assemble an arbitrarily large number of qubits in a cat state in a time much shorter than T2 or, 
more easily, by assuming that one can take as long as desired to determine g using an arbitrarily 
large number of qubits prepared in a product state. Once decoherence becomes a consideration, 
the uncertainty in estimating g should be written in terms of the decoherence time T2 and the 
relevant resources, R and r, not directly in terms of the number of qubits used. 

The results of our analysis, summarized in Table [J and Fig. [3l show that cat-state entanglement 
is useless if r > T2. When T2 is much larger than r, the results show that one should put as many 
qubits as possible in each cat-state probe. For intermediate decoherence, it is not that cat-state 
entanglement is useless, but rather that one should make a judicious, optimal choice of how many 
qubits to include in each cat-state probe. The overall conclusion is that entanglement is only 
useful when one can make the effects of decoherence small on the time scale over which one must 
estimate g. While this conclusion is reached here for a special model of measurements on qubits, 
it is generally true for quantum-limited measurements in the face of decoherence. 

Our analysis highlights one further point, having to do with using the quantum Cramer-Rao 
bound to determine quantum limits. As long as one is interested only in measurements involving 
pure states, the form of the Cramer-Rao bound as a generalized uncertainty principle is sufficient 
for investigating bounds on measurement accuracy. Once decoherence is introduced, however, 
inevitably leading to measurements on mixed states, one must use the stronger form of the Cramer- 
Rao bound involving the Fisher information to obtain meaningful bounds on measurement accuracy. 
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APPENDIX A: GENERALIZED UNCERTAINTY RELATIONS FOR PRODUCT STATES 

Consider a continuous trajectory in the space of product states of n qubits parametrized by X, 

n 

p p (X)=® Pj (X), (Al) 
i=i 

If the value of X changes by a small amount dX, the change in pj can be written as 

ihidX „ AhjdX 



pj -> pj + dXp'j = J2 dp a] \aj){aj \ + e ih ' dX pje 

The vectors make up the eigenbasis of pj, with eigenvalues p aj , i.e. 

Pi = £p«il«i><«il • 



(A2) 



(A3) 



The operators hj are the generators of translations in X for each of the n systems, while dp aj are 
small changes in the eigenvalues of pj due to the small change in X. Notice that in the presentation 
of Sec. HH the fact that At is independent of the parameter g means that the eigenvalues of p(g, T) 
do not change with g; thus the terms having to do with eigenvalue changes do not appear in that 
discussion. 

Keeping terms to linear order in dX, we have 



(A4) 



The corresponding change in the overall state, 



P P ^ P P + dXp' p = p p + dXY^ Pj Pk 

j=i k^j 



(A5) 



gives us 



p' P = 12 

3=1 



)Pk ■ 



(A6) 



Our objective is to obtain an expression for a line element dSp in the space of density operators 
p p that measures the distinguishability of neighboring quantum states. Following Eq. (12. 5|) . we 
have 

2 



dspY 
dX 



tr 



Pp£p P (p'p) 



(A7) 
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We start by computing C Pp (pL) using the definition in Eq. f|2.6|) . Noting that p p is diagonal in 
the tensor product basis furnished by {|ay)} for each of the systems and using Eq. (|A6j) we obtain 

2 

CPP{PP) ' ^..^5]?..^ PPl-'-PPn+PSl-'-PSn 



x(p 1 ...p n 



E 



Mi 



«5i...«5n^|A.../^)(5 1 ...«^| • 

(A8) 



Simplifying this expression gives 



C pM'p) = E 



E (dPaj/dX) i^ 3 -i + 2«E Paj " Pft 



Pc 



Pay+Pft 



Mi 



Mi 



(A9) 



which leads to 



a=l 



i|) = E tr[^« (#] x II <*P* = E tr[^„ (#] = E (if) • ( Al 0) 

j=l fc^j j=l j=l 

where we use tip'j = 0. In the special case where all pj are identical and equal to p and when 
changes in X affect all the systems in the same way, we have 

2 



dsp f— ds . ( ds \ 

— — = vti where = tr 

dX V dX \dX 



(All) 



APPENDIX B: USEFUL IDENTITIES INVOLVING PAULI OPERATORS 

A few identities involving products of Pauli operators that are used in our calculations are listed 
below: 



J2°z;j® t k)®(t j ±a z;j ) = ±n(g){l j ±a 

3=1 



z;jj j 



\i =1 J+k J j=l 
n \ n 

y=l j^k J j=l 



i a 



E a *-j ®°z;k h 



\3+k=l 



W,k I 3=1 



3=1 

n 

in 2 - n) (g) (% ± a z . tj ) 

3=1 



^2 a z -j ® a z . k (g) lA (a x -j ± ia y .j) = (n 2 - n) (g) (a x .j ± icr^y) , (Bl) 

\j^k=l l^j,k J j=l 3=1 



E ® Ife C8) (% ± Avzu) = ±^E ± <S> C 1 * ± ^i*) > ( B2 ) 
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®(1,±A7 Z;J ). (B3) 
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